ON THE BLOWING UP OF SOLUTIONS TO QUANTUM 
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Abstract. The blow-up in finite time for the solutions to the initial-boundary value problem 
associated to the multi-dimensional quantum hydrodynamic model in a bounded domain is proved. 
The model consists on conservation of mass equation and a momentum balance equation equivalent to 
a compressible Euler equations corrected by a dispersion term of the third order in the momentum 
balance. The proof is based on a-priori estimates for the energy functional for a new observable 
constructed with an auxiliary function, and it is shown that, under suitable boundary conditions and 
assumptions on the initial data, the solution blows up after a finite time. 

Keywords: Quantum hydrodynamic equations, blow up of smooth solutions. 

1. Introduction. The evolution of a quantum fluid in a first approximation can 
be described by a dispersive perturbation associated to the Hamilton-Jacobi system 
for compressible fluid dynamics, sometimes referred as a dispersive perturbation of 
the Eikonal equation for the evolution of amplitude and phase velocity of quantum 
wave guides. This system consists in the continuity equation for the particle density 
p and for the momentum pu 

p t + div(pu) = 0, OO.ieSJCl' 1 , (1.1) 

(pu) t + dW(pu ®u) + VP(p) = » (1-2) 

subject to the initial conditions 

p(-,0) = pj(-) >0, u(.,0) =«/(•), (1-3) 

where P(p) > describes a pressure function, due to the boundary or to mean field 
effects, and e denotes the scaled Planck constant. 

Formally equations (jl.ip . (|1.2|) can be derived through the Madelung's transform [12] 
from the nonlinear Schrodinger equation 

^t = -yAV> + MM 2 )^, xeR d , (1.4) 

i>(;0) = M-) (1-5) 
for the wave function ip(x,t), where h is an integrable function such that 

P 
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Equation (|1.4p has been proposed as multi-particle approximations in the mean field 
theory of quantum mechanics, when one considers a large number of quantum parti- 
cles acting in unison and takes into account only a finite number of particle-particle 
interactions. In fact the complex- valued wave function ip can be rewritten as 

i\> = Vpexp 

where p = \ip\ 2 and S is some phase function. Then by plugging the above formula 
for the wave function to (|1.4|) the imaginary and real parts of the resulting equation 
give respectively (|1.1[) and (|1.2|) with u — VS. The formal equivalence between (jl.ip . 
(|1 .2|) and (jTT4j) is explained more in details in the following lemma: 

Lemma 1.1. Let if) be a solution to hi .5]) with initial datum ipj := ^/pjexp 

the functions p := \ip\ 2 and u :— e Im ^^^ solve U.1\) - H77S\) with initial data pj — \ipi\ 2 
and ui = VSj as long as p > 0. 

Vice versa let (p, u) be a solution to ll-l]) - i TO)) . where u = VS" and p > in Q, t > 0; 
the function ip := y/pexp (^) solves {1.4]} • 

In the recent years problem (|1. 1|) . (|1.2p has attracted a lot of interest and many 
papers have been published. In [lOj Jiingel and Li proved local-in-time existence 
of smooth solutions to (jl.ip , (| 1 . 2[) for the one dimensional case with Dirichlct and 
Neumann boundary conditions for the particle density p. Moreover, if the initial data 
are close enough to the steady-state solution, the local-in-time solutions are shown to 
exist globally in time (so called small solutions) . 

In [3] j 0] the authors studied the positive steady-state solutions to the potential flow of 
(|1.1[) , (jl.2p . In that case the one dimensional stationary problem can be reformulated 
into the following form 

pu=J, \u 2 + h{p)-^=K, (1.7) 

with appropriate boundary conditions, where J and K are given constants (in the 
classical theory of compressible hydrodynamics, J is the mass flux, K the isoener- 
getic constant and the function h denotes the enthalpy function). For this Quantum 
Hydrodynamic case the authors proved in [3], [4] that, if K is sufficiently small and 
J large enough, system (|1.7p can not have even a weak solution (supercritical con- 
dition). In addition they showed that one can not recover the entropic solutions of 
the compressible Euler limit for the supercritical case unless the model incorporates 
diffusion effects. This suggests that a smooth solution to the corresponding evolution 
problem (|l.ip . (|1.2p . for which the iso-energy functional K, defined as 

K(x,t)= Q" 2 + Mp)-^|f) (x,t), (1.8) 

is small at the boundary of the domain, will blow up in finite time. This result is 
proven in Theorem 12.11 in Section [3 and generalized for a multi-dimensional domain 
in Theorem 13. 11 

The existence of a one-dimensional stationary viscous quantum hydrodynamic model 
is shown in [8]. This model corresponds to a quantum regularization in the Fokker- 
Planck collision operator at the Wigner level [1] , which generates a linear second order 
viscous term in (jl.lj) and (|1.2p , after a moment method [S] . This regularization allows 
for solutions in case when a "weakly" supercritical condition holds. The existence of 



solutions in the supercritical regimes for this model is still an open question. 
In [5] the author proved non-existence of non-constant traveling waves in the super- 
sonic regime for the Gross-Pitaevskii equation (|1.4[) with <7'(|-0| 2 ) = (1 — \ip\ 2 ) in 
any dimension d > 2; in particular it is shown that any quasi-stationary solution 
in the whole space of the form ip = ip(x — ct,X2, ■■•,%d)i with finite energy and ve- 
locity c > V2, is constant. The proof is based on integral identities for the variable 
rj := 1 — |^| 2 and estimates for the energy functional E(ip) := J Rd | VV>| 2 + (1 — IV'I 2 ) dx 
using the Pohozaev identities. 

Local in time existence of solutions of (|1.2[) in the whole space has been shown 



Our method of showing blow-up for a boundary value problem associated to the 
system (jl.ll) . (II. 2p is in essence inspired in the original argument of R. Glassey [TJ, 
where he proved finite time blow up of smooth solutions to the focusing nonlinear 
Schrodinger equation (| 1 .4[) with large initial data. The main idea of the proof in [7J is 
the study of the nonnegative quantity I(t) = J Rd \x\ 2 \tfj(x, t)\ 2 dx. More precisely the 
author showed that lit) will be negative if the smooth solution to (|1.4p exists long 
enough for the focusing case. The observable I[t) becomes close to zero due to mass 
concentration at the point x = 0. 

This idea has also been used by Sideris [H] in the study of the finite time blowing up 
of smooth solutions to the compressible Euler equations and by Xin [15j in the study 
of compressible Navier-Stokes equations. 

In our work the blow up of the system is caused from generation of vacuum almost 
everywhere due to boundary constrains. One may expect (or conjecture) this effect 
is caused by mass concentration at the boundary of the fluid domain. Such study of 
the type of singularity formation remains an open question and it is not addressed in 
this manuscript. 

Since we are dealing with a boundary value problem, the weight function a(x), which 
defines the new observable I(t) = J n a(x)p(x, t) dx, should be chosen differently from 
the Cauchy-problem case. It will turn out that the function a(x) must be a concave 
function with zero value on a part of the boundary and its form will depend on the 
domain and on the boundary conditions of the problem. 

This paper is organized as follows: in Section[2]we investigate the finite time blow- 
up of smooth solutions to the one-dimensional quantum hydrodynamic equations (we 
assume for simplicity e 2 = 2) 



in [TT]. 



Pt + (pu) x = 0, 



t>0,i£ [0,1], 



(1.9) 




(1 



10) 



with the initial and boundary conditions 



p(;0)= Pl (-)>0, U(;Q) 
Px (0,t)= Px (l,t)=0, 



(1 
(1 



11) 

12) 




(1 



13) 
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where C\ and c 2 are negative constants. Assumption (|1.13p with c\, c 2 < means 
that there is a nonnegative total momentum flux across the boundary, which implies 
small values for the iso-energy functional K(x,t) at the boundary for all time. 
Under Neumann () 1 . 1 2[) and Dirichlet boundary conditions for the particle density p 

p(0,t) = P!, p(l,t)=p2,, pi, P2 > 0, 

local in time existence of classical solutions and global existence of small solutions to 
(|Q ]) -([UT 1) is proved in [TU]. 

We point out that the proof of local existence of smooth solutions to (|1.9[) - (|1.13[) re- 
mains at this stage of the work an open problem and will be addressed in future work. 
We believe however that the same techniques used in [10] can be applied for this case 
and existence results will read as follows: if the mass flux J = pu at the boundary is 
small enough and the momentum flux across the boundary is large enough (ci , ci > 
in (|1.13[) ). system (|1.9p - (|1.13p will have global in time existence of solutions. On the 
other hand the authors in [3] proved that if the total momentum flux is small enough 
and the mass flux is large, (ci = c 2 < 0) problem (11.9| - (|1.13|) does not have stationary 
solutions. Our conjecture therefore is that the proof of existence of smooth solutions 
to (|1.9|) - (|1.13|) with a, C2 negative will be possible only locally in time. The results 
shown in this paper address the fact that, in the latter case, any smooth solution 
to (|1.9p - (|1.13p will blow-up in finite time. However our results do not exclude that 
weaker solutions of the problem may exist for longer time; our interest is related only 
on solutions (p, it) that deliver strictly positive particle density everywhere in the fluid 
domain. 

Following the same idea mentioned before for the whole space problem in the Schrodingcr 
approach, we are able to show that in a bounded domain problem (|1.9|) - (|1.11|) sub- 
ject to the boundary conditions (|1.12p . (| 1 . 1 3[) does not have global- in-time classical 
solutions: more precisely it will be proved that under certain assumptions on the 
initial data, the nonnegative quantity I(t) = a{x)p{x,t) dx with a(x) > 0, for 
all x E [0, 1], will become strictly negative after a certain finite time T* , where T* 
depends only on the initial and boundary data of the problem. That implies that 
smooth solutions p can not exist any longer for t > T*. In this sense we interpret 
finite time blow-up of a smooth solution. 

The proof is based on two important points: the first one relies on a-priori estimates 
for the total energy functional 



consisting on kinetic energy, thermodynamical total enthalpy and quantum energy of 
the system, where g(p) :— J P h(r) dr. The second point is related to some assump- 
tions on the iso-energy functional K(x,t), defined as in (|1.8p . which is bounded by a 
nonpositive constant at the boundary of the domain for smooth, local in time positive 
solutions p and u. Note that (|1.8p becomes the classical Bernoulli's law when the 
dispersive term is absent. It is our aim to stress that the role of the Bohm potential 
in K is crucial in order to prove lack of existence of solutions globally in time. In the 
classical fluiddynamical framework, the third order dispersive term derived from the 
Bohm potential is absent in the conservation of momentum equation and it is well 
known that in this case the system has well-posed solutions for any time. 
Section [3] shows how to extend the same method to the multi-dimensional equations. 
In the last section we give an alternative configuration for the one-dimensional case 




(1.14) 
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where system (|1.9p - (|1.10l) does not support smooth solution for any time. 
We show how the weight function a(x) used in the definition of the observable 1(f) 
depends strongly on the domain Q. The computations are shown in a general domain 
and two examples are giving at the end of the paper. 

2. Blow-up in one-dimensional space. In order to present in a simple manner 
the idea involved, we first consider the problem in a one dimensional domain given 
by a simple interval. The main results of this section are formulated in the following 
theorems: 

Theorem 2.1. Let p € # 1 (0,T,.L 2 (0, 1)) n L 2 {0, T, H 3 {0, 1)), p > in [0, 1] x 
[0, T], lie H 1 ^, T, L 2 (0, 1)) n L 2 (0, T, H 2 (0, 1)), be a small solution to fOP - ITX^ 
with initial conditions < pi € -ff 1 (0, 1), uj 6 L 2 (0, 1). 
Let uj and pj be such that 

»i 

(1 - 2x)p I u I dx = M < 0. (2.1) 



There exists a positive constant t$ < T* < +oo depending only on the initial and 
boundary data such that for t > T* any classical solution p does not exist any longer. 
In particular lim^r* p(x, t) = almost everywhere. 

Remark 2.2. One may interpret in the above theorem that classical solutions p 
lose regularity when t approaches T* . This is our interpretation of finite time blow-up. 
Theorem 12.11 is expanded to higher dimensional sets in the next section, also under 
some existence conditions on the auxiliary function a(x). 

In order to prove the above theorem, we need some preliminary results. 
Lemma 2.3. Let u, p be smooth solutions to hl.9\ )- Fl.l3\) and E(t) and K(x,t) 
be defined as in {1. 1J$ and hl.8\) respectively. For every t > it holds 

E(t)-E(0)+ [ (upK)(l,s)-(upK)(0,s)ds = 0. (2.2) 



Proof. Taking <fi — u as test function in (|1.10| and integrating with respect to x 
over [0, 1] we get 



1 fl r 



(up) t u dx — 



o 



u(pu 2 ) x + u(P(p)) x - up 



( Vp xx 
V VP 



dx. (2.3) 



We multiply now (II. 9|) by the test function <f> = h{p) — — \u 2 : 

\t[h{p)-^- l -u 2 ") <lr (2.4) 







/ [up) x ( h(p) - - -u' ) dx. 



Since p x = on the boundary, it holds 



J [up) t u - ^p t u 2 dx = J ^pu 2 dx, 
J h(p)p t dx = ^ t J 9{P) dx, 

1 Vp xx f\ ^ 2 



Pt ^ dx= diJ ^ dx - 
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Therefore, after adding up the left- and the right-hand side of (|2.3p and (|2.4p it follows 
9E ^ - I u(pu 2 ) x dx- f u(P(p)) x dx+ f upf^A dx 



or Jo Jo Jo \ y/ 5 /x 



(up) x h(p)dx + / (ttyo)-,, — + - / (li^u da; 

VP 2 Jo 



= :/l+/ 2 +/ 3 +/4+/5+/6- 

Noticing that 



i r 1 

■fi + 4 = -g y (/ou 3 )^ da:, 

h+h = lo\ PU ^ff) 
h + h = - \ {pu h(p)) x dx, 



the thesis follows. 



□ 



Lemma 2.4. Letu, p be smooth solutions to (1.9\) - (T.l$fy andl(t) be a nonnegative 
quantity defined as 

I(t) := [ x(l-x)pdx, Vt>0. 
Jo 

It holds 

il{t)= I (1 - 2x)piui dx - 2 / f {pu 2 + P(p) + 2y/f x ) dxds 
at Jo Jo Jo 



+ p lu 2 + x -^-^\(l,s)ds 



p Vp 

+ fpU+^-^A (o, s) d S . 

Jo V P VP / 

Proof. Integrating (|1.10|) with respect to x and i over [0, 1] x [0, t] we get 

pu dx — / piuj dx + p ( u 2 H I (1, s) ds (2-5) 



f p ^ + P(£)\ (0j fi) ds = (^/p xx (l, s ) - ^p^p xx (0, s)) ds. 

We multiply now (jl.lOp by a; and integrate over space and time; taking into account 
that 

J J xf} (~^j dxds = J ^/p^/p xx (l,s) ds - J J xp x ^^- dxds 

pt pi pt pt pi 

VPVP XX dxds = / VPVPxxi 1 ' s)ds + 2 / yfp 2 x dxds, 
Jo Jo Jo Jo Jo 



it holds that 



/ xpudx— I xpiui dx — I I (pu 2 + P(p) + 2s/p x ) dxds (2-6) 
Jo Jo Jo Jo 

Now we take (j> = x(l — x) as test function for (|1.9[) : 

/ x(l — x)pt dx — I pu dx — 2 I xpu dx. (2-7) 
Jo Jo Jo 



d_ f 1 

di. 



Substituting and in we get 

/ x(l — x)p dx = / (1 — 2x)pjuj dx — 2 / / (pu 2 + P{p) + 2^/p 2 ^) dxds 
Jo Jo Jo Jo 

+ rp(u 2 + ^i-^)(i, s) ds 

Jo I P Vp J 

and the thesis follows. □ 

Proof of Theorem \2.1\ Directly from Lemma [2T4l using (| 1 . 1 3[) and (|2.ip . we get 
I(t) < 7(0) + 1 f (1 - 2x)uipi dx = 1(0) + M t, 



where 1(0) — J x(l ~ x)pi dx. It is easy to see that I(t) < if t > T* :— 
which implies blow-up for the solution p at time t = T*. □ 



3. Multi-dimensional generalization. In this section we extend some aspects 
of the method developed in the previous section for one space dimension to any space 
dimension, provided the boundary value problem has a solvable auxiliary problem in 
the fluid domain Slcl^M. 

Our argument consists in constructing a nonnegative weight function a(x) defined 
over the set Ct such that 

Hess(a) is negative-semidefinite in Q, (3.1) 

Aa = -g(x) < in Q, (3.2) 

a = on dfl D C dtt, (3.3) 

Va • v = on dfl N = dfl \ dVL D , (3.4) 

for which the new observable 

I(t) = I a(x)p(x,t) dx 
Jo 



becomes negative in finite time. The function g(x) must be nonnegative and satisfy 
the following conditions 

Ac < in a ^ > on dVt. (3.5) 

Examples of domains where problem ()3. (|3.5j) can be solved are given in Remark [3.3 1 

In the following theorem we show that the auxiliary weight function a{x) satisfying 
(|3.ip - (|3.5jl is good enough to create an observable quantity for which the d-dimensional 
quantum hydrodynamic system (|l.l[l - (jl.3|) will admit solutions that blow up in finite 
time. 

Due to the constrains on the existence of the auxiliary weight function a(x), the 
system (|1.1[) - (|1.3I) under consideration is subject to the following boundary conditions 
posed on the fluid domain f2: 



da /A^p | |Vyp| 2 P(p) 
dv\^fp p p 



-(u-Va)(u-u) = c<0, ondn D , (3.6) 
u-v — 0, ondCl N , (3.7) 

dv 



0, on dfl, (3.8) 



where dtt = dflo^dflN, v is the outward normal unit vector of dfl and c a nonpositive 
constant. The boundary condition Q3.6P can appear to the reader in a first moment 
not completely correct, since the weight function a(x) appears in its formulation and 
depends on the domain. We will show in Remark 13.31 that for sufficiently smooth 
domains, the geometry of a(x) is such that condition (|3.6|) contains exactly the same 
information as boundary condition (|1.13|) for the one dimensional problem (which 
means nonnegative flux across the Dirichlet boundary 8Qd)- 

Existence of solutions to (|l.ip - (|1.3p is still an open problem. The difficulty in the 
existence proof is given by the complexity due the multidimcnsionality of the problem. 
However non-existence of stationary solution to (|1.1|) - (11.3|) . (|3. 6|) - ()3. 8|) is shown in [3]. 



Theorem 3.1. Let p e H 1 (Q, T, L 2 (fl)) n L 2 (0, T, H s+2 (n)), p > in Q x [0, T], 
u £ ff 1 (0,r,I 2 (n)) n L 2 (0 : T,H s+1 (fl)), for s > f, be a solutions to rO)-fOI). 
i3.G\) - nF8\) . Leta(x) > be a nonnegative function satisfying hS. J)) -i TO|) . If pi andui 
are such that J n piu\ ■ Va dx < 0, there exists a finite time to < T* < +oo depending 
only on the initial and boundary data such that for t > T* any classical solution p 
does not exist any longer. In particular lim t ^T* p(%, t) = almost everywhere. 

Remark 3.2. This is the extension of Theorem \2.1\ to higher dimensions. 

Proof. We consider the quantity I(t) = j n a(x)p(x,t) dx, where a satisfies (|3.1[) - 
(|3.4|) , (|3.5p . Using (|1.9p and taking into account that a = on 8Qd and u ■ v = on 
drift, the time derivative of I(t) reads as follows 



9I(t) 
dt 



/ ap t dx = — / a d\v(pu) dx 
Jn Jn 

apu ■ v ds + / pu ■ Va dx 
an Jn 

pu ■ Va dx, (3-9) 

o 
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where v is the outward normal unit vector of dfl at x. Let Va be a test function for 

/ (pu)t • Va dx + Va • div(pit ® u) da; + / Va ■ VP(p) dx 
Jo Jo Jo 

- ^ pV ■ Va dx = =: /i + J 2 + 7 3 + 4- (3.10) 

Since s > § then p(t, •) € H s+2 {n) ^ C 2 (ft) and «(t, •) € H s+1 (n) ^ C 1 ^) arc the 
classical Sobolev cmbcddings with the L°°-norm in the target space. Integrating by 
parts it holds 



h = Va ■ div(pu eg) u) dx 
Jn 

— I div(pu ■ Vau) dx — pu T Hess(a)u dx 
Jn Jn 

= / p(u ■ V a) (it ■ v) ds — / pu T Hess(a)u dx, 

Jon D Jo v v ' 

<o 

taking into account (|3.7p . In addition, after integration by parts, we get 

h= [ Va • VP(p) dx = [ P{p)^- ds - [ P(p)Aa dx, 
Jn Jdn ° v Jn 

and 

U = - ( pV f — ^ ] • Va dx 



sTP J 

[ ^/pAy/p^- ds + 2 [ A^/pVy/p ■ Va dx + [ yfpAyfpAa dx 
Jon ov J a J n 

[ VP A %/P7T ds - 2 I V VP ' v ( Va ' v Vp) dx + I { 9\ V V?| 2 + ■ v 5 J da; 
Jan 0t/ Jo Jo V 1 J 

/ ^fpA^fp^- ds-2 V y/p{Ress(y/p)V a + Hess(a)V v /p) dx 
Jeo Jn 

[ 9\^\fp\ 2 dx + \ f p^- da- \ ( pAg dx 
Jn z Jon ov i Jn 

[ ^/pAy/p^- ds - [ V|V^/p| 2 ■ Va dx - 2 f V ^/p T iless{a)V ^/p dx 
Jon ov Jn Jn 

[ S|VV?| 2 dx + i / P^ds- \ ( pAg dx 
Jo 2 Jao 2 Jo 

Jen ^ A ^° + I V ^! 2 ) - \ p ^£) ds - 2 J n V\/p T Hess(a)V^ + ^pAg dx, 



<o 

using (|3.2p , (|3.5p and (|3.8p . Substituting the above computations in (|3.10p we obtain 

/ {pu)t ■ Va dx < pu Hess(a)w dx+ / (P(p) + 2|VVp| 2 )Ao dec 
Jo Jo Jo 

-da (A^p |VVp| 2 P(p)^ 



+ / P 

ldn D 



dv V Jf> 
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— (u • Va)(u • z/) 



ds. 



This implies 



pu ■ Va dx < piui ■ Va dx. 
n J Q 

Integration of (|3.9p with respect to time and the above inequality yield 



I(t) = 1(0) + 



/ pu ■ Va dxdr < 1(0) +tl piui ■ Va dx. 
Jo Jn Jn 



If pi and uj are such that J^piui ■ Va dx < c < 0, then for t > — j pj ^°y a it 
holds J(i) < and the thesis follows. 

□ 

Remark 3.3. We finish by showing examples of domains and observable I(t) in 
which, if 13. 6\) holds, any solution to i i.lj )-i TOj) . 13. 6\) - %3J$j) blows up in finite time. 
The construction of the weight function a(x) satisfying H3.1]) - [3~5\) depends on the 
shape (geometry) of the fluid domain and on the boundary conditions of the system of 
equations associated to it. 

Constraint \3.5}) is satisfied in particular when g = c, with c a nonnegative constant. 
In this case a solution to A3. l\) - (3~^\ ) can be given by a quadratic polynomial function 
or by a solution to the Monge-Ampere equation 

det(Hess(—a)) = 1, a = on <9f2, 

on a smooth convex domain, such that Hess(a) is a diagonal matrix (see £|J/ or 
for a general survey in the area). 

We give two easy examples of domains. 

First case: spherical domain. We consider a spherical domain := {x £ R d | \x\ < 
1} for problem (|l.lj) - (|1.3|) with boundary conditions (|3.6D - (|3.8[) . In this case the 
function a is given by a(x) = 1 — \x\ 2 , x £ 57. It holds u T Hess(a)u = — 2|u| 2 , 
Aa = —2d on f2 and |^ = _2 on dfl. Condition (|3.6| reads as follows 



P 

on 



da /A> P(p) 
dv\^fp p 



(u ■ Va)(?i • v) 



ds 



2 / P 

Jon 



(u ■ x) z + 



2 , P(p) avp' 



\[P 



ds, x S dfl. 



Remark 3.4. The one dimensional case corresponds to this example if the inter- 
val considered is [— 1, 1]. 

Second case: cylinder-like domain. Let Q be a bounded cylinder-like domain 
R d D n := [-1, 1] x Q.1, where D,i C M^ 1 . We denote with z = (x,y) for x G [-1, 1] 
and y G fli- 

In this case the boundary dflo = fli X {x = —1} U ill x {x = 1} and dflN — 
90 x x(-l,l). 

We choose a(x,y) = 1 — x 2 ; it holds Aa = —2 and u T Hess(a)u = — 2u 2 on f2, where 
u x denotes the component of the velocity along the x-direction. In this case (|3.6p 
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becomes 



P 



ill x{x=— 1} 

P 1 K 



p(pY 


) -(«■ 


p / 


p(p) 




p 


/ 


p(p) 




p 





P i u l + — — - — 7=- ) (0,y,s)dyds 



ds 



f!ix{x=l} 

taking into account that w • 1/ = and Va ■ 1/ = on <9f2 at . 

Remark 3.5. The extension of Theorem \4- l\ to higher dimensional spaces requires 
more delicate assumptions in the existence for the initial boundary value problem and 
corresponding extensions to d-dimensions for ^4.2\ ). It is expected that the same kind 
of result will follow for the two cases of fluid domains covered in this section. 

4. Another configuration for possible blow-up. In this section we give an- 
other possible blow-up configuration of solutions to (| 1 . 9[) - (| 1 . 10[) under different con- 
ditions at the boundary. We consider the system in one dimensional space subject to 
zero Neumann boundary conditions for the particle density (|1.12|) and Dirichlct-type 
conditions for the velocity 

u(0,t)=u o , u(l,t)=u 1 , Vt>0. (4.1) 

Local in time existence of solutions of such problem (|1.9|) - (|1.10|) . (I1.12| . (|4.1|) can 
be shown using the same technique as in [lOj . The problem of global in time well- 
posedness for the system is related to the lack of maximum principle for the current 
density J = pu. Formations of vacuum for the particle density implies no bounds 
for the velocity of the fluid. We believe that smooth solutions to (|1.9j) - (|1.10p . I|1.12[) . 
(|4.ip with initial data such that -E^O) is large enough and K(0) is sufficiently small 
(where E and K are defined as in (|1.14jl and (|1.8p respectively), will deliver also in 
a later time small value for K(-,t) at the boundary of the domain (condition (|4.2p ). 
The rigorous proof of such statement is still an open problem. 

Theorem 4.1. Let p € fT 1 (0,T,.L 2 (0, 1)) n L 2 (0, T, H 3 (0, 1)), p > in [0, 1] x 
[0,T], ue ^(0,7,^(0, 1)) n L 2 (0,T,H 2 (0, 1)), be a solution to fOP- flJl)) . glp 
where P(p) satisfies |^.g[ ) and &4-4\ l- Let < pi E H 1 ^, 1), Uj £ L 2 (0, 1) be such that 
E(0) = Jp 1 \piu 2 + g(pi) + (\/pT)x dx is large enough. If the quantity K(x 1 1), defined 
as in U.8\) . at the boundary of the domain is negative and bounded in the following 
way 

-M < K(0,t) < -amax 2 {\u \, |tti|), for all t S [0, T], (4.2) 

for some constant M > and a > 1, then there exists a positive constant to < 
T* < +oo depending only on the initial and boundary data such that, for t > T* , any 
classical solution p does not exist any longer. In particular limt^T* p(x,t) = almost 
everywhere. 

For the next theorem we need some assumptions on the pressure function P(p) 
of the system. More precisely we assume that 

P ^ > A > 0, V p > 0, (4.3) 



g(p) 
p(p) 



Hp) < o, v p > o. (4.4) 



n 



Remark 4.2. Assumptions {-j-Sty and f^.^p are satisfied for example by P(p) = p 1 
for 7 > 1 and by any other sum of power functions. In this case g(p) — t^ziP 1 an d 

9(P) - 1 

Proof of Theorem\Jl\ We first show that K(0, t) = K(l, t) for all t > 0. This is 
expected since K(x,t) becomes the "isoenergetic" constant if stationary states hold 
and it is a simple consequence of the conservative character of the underlying Euler- 
type system. In fact, without loss of generality, we can assume that 

u(x,t) = S x (x,t), S(0,t) = C o , S(l,t) = d. 

This is the classical assumption to writing down the quantum hydrodynamic system 
from the Madelung transforms applied to the wave function, solution to the corre- 
sponding Schrodinger equations (see Lemma ll.ip . 

Therefore f Q u(x, t) dx = C\ — Co- Dividing now (I1.10|) by p and using ljl.9|) we get 



2 P v J, 

Integrating the above equation with respect to x over [0, 1], it holds 

d f 1 

= —/ u(x,t) dx = K(0,t) -K(l,t), 
ot Jo 

which implies K(0, t) = K(l, t). 

Next, taking into account (14. 4j) . the following inequality holds 



.2 , P(P) VP XX \ s J\.2 



p[u 2 + -^-^)<p^ + K(x,t) 



and therefore from Lemma \2. 41 we get 



I(t) < I (1 - 2x)p I u I dx - 2 / [ (pu 2 + P(p) + 2y/pl) dxds 
Jo Jo Jo 

- f p(1,s)(k(1,s) + 1-uI) ds+ f P (0,s)(k(Q,s) + ^u 2 ) ds, 



d_ 

di 



where I{t) = f Q x(l — x)p dx. We assume that max(|uo|, |ui|) 7^ 0. 
Using (14. 2|) on the above inequality we get 

^I(t) <M Q -2jJ^ (pu 2 + P(p) + 2jp 2 x ) dxds 

a - max 2 (|u |, J p(0, s) + p(l, s) ds, 

where Mo '■= J (1 — 2x)piui dx. 

If Mo < the thesis follows as in the previous theorem. 
If Mo > we divide the proof in two parts: 
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(4.5) 



First part: If either f*p(0,s) ds > 2M ° , or f*p(l,s) ds > 2M °, , for 

t > 0, it holds 

I{t) <I(0)-M t(2a-2), 
where 1(0) = L x(l — x)pi dx, and the theorem is proved. 

Second part: If both J* p(0, s) ds < and J*p(l,e) ds < ^f^M) ■ 

then we need the estimates from Lemma 12.31 and it follows 



E{t)=E(0)+ f K(0,s)(p(0,s)u o - p(l,s) Ul ) ds 
Jo 

>£(0)-max(h)|,kl) / \K(0,s)\(p(0,s)+p(l,s))ds. 
Jo 

We recall that Lemma 12.31 can be proved also for solutions to l|1.9|) - (|1.10|) . (|1.12p . 
(JUT]). Using now to control K(0,t) we get 

AMnM 

E(t) > E(0 . 

max(|u |, |«i|) 

Assumption (14.411 implies that 

-2 / / pit 2 + P(p) + 2(^) 2 efecZs < -2min(A, 2) / £(s) ds , 
Jo Jo Jo 



and finally from the differential inequality (|4.5[) . it follows 

4M M 



J(i) < 7(0) + M t - min(A, 2)i 2 E(0) 



max(|u |, |ui|) 



Therefore if E(Q) > m J^\ Ull) and t > T* , where T* depends only on the initial 

and boundary conditions, we get x(l — x)p dx < for all t > T* and the theorem 
is proved. 

If uo = u\ = 0, then 

x{l - x)p dx < [ x(l - x)pi dx + 1 (M - 2£(0)min(A, 2)) , 
Jo 

and therefore if E is large enough such that E > „ M ( a . „. the thesis follows once 



+ \ T^* Io 

L ^ ± -~ 2B mm(A,2)-M " 



2min(A,2) 

□ 
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